icin sina LAN KAA ULLULhdntainadnnnunnnARLEn sinuses ensue 
q 
D 


DISCLAIMER NOTICE 


THIS DOCUMENT IS BEST QUALITY 
PRACTICABLE. THE COPY FURNISHED 
TO DTIC CONTAINED A SIGNIFICANT 
NUMBER OF PAGES WHICH DO NOT 
REPRODUCE LEGIBLY. 


Ap? 27001 


AL 71-006} 
MARCH 197% 


Aerespace Research Lahoratories 


es —_——— 


THE DESIGN OF AN AXIAL COMPRESSOR STAGE 
FOR A TOTAL PRESSURE RATIO OF 3 TO 1 


3. I. WENNERSTROM 
RICHARD 14. HEARSEY 


FLUID DYNAMICS FACILITIES RESEARCH LABORATORY 


PROJECT NO. 7655 


Approved for public release; distribution unlimited. 


= AUG 2 197% 
AIR FORCE SYSTEMS COMMAND f he 1 
United States Air Feree (0507, 6 
Resredred By cs 
NATIONAL TECHNICAL xz 
INFORMATION SERVICE 
Seen-fed va 2181 


\Zu 


NOTICES 


When Government drawings, specifications, or other data are used for any purpose other than m 
‘connection with a definitely related Government procurement operation, the United States Government 

‘ thereby incurs no résponsibility nor any obligation whatsoever; and the fact that the Government may 
have formulated, furnished, or in any way supplied the said drawings, specifications, or other data, is 

_ not to be regarded by implication or otherwise as in any manner licensing the holder or any other 
person or corporation, or conveying any rights or permission to manufacture, use, or sel] any patented : 
invention that may in any way be related thereto. : 


We ANAND Ea URI a at he IE ke ew 


Bag ae IN Ors 


Agencies of the Department of Defense, qualified contractors and other 
- government agencies may obtain copies from the 


Rar pe Cera ner coears 


Defense Documentation Center 
Cameron. Station 
Alexandria, Virginia 22314 


seks Ft 


This document has been released to the 


Sue LENA ene? be 


CLEARINGHOUSE | 
payouts o U.S. Departrnent of Commerce 
q Springfield, Virginia 22151 


ISTIIEATION narsennnen at 
for sale tojthe public, 


ea cunnenentansscesesoconsmsoneeonaponeoue eet 


’ Copies of, ARL Technical Documentary Reports should not be returned to Acrospace Research 


Laboratories unless return is required by security considerations, contractual obligations or notices on 
p a specified document. 


200 - June 1971 + CO305 - 42~71-723 


Re atta, tee 


% 
iy 
EF + 
Re 
me 

43 


es 


UNCLASSIFIED 


security Classification 


DOCUMENT CONTROL DATA-R&D 
(Security classilication of title. body of abstract and indexing annotation must be entered when the overetl repost is classified: 
Ri ORIGINATING SCTiviTty (Corporate author) 


Aerospace Research Laboratories (ARL/LF) 
Wright-Patterson AF Base, Ohio 45433 


l2e. REPORT SECURITY CLASSIFICATION 


Unclassified 


US. REPORT TITLE 


THE DESIGN OF AN AXIAL COMPRESSOR STAGE 
FOR A TOTAL PRESSURE RATIO OF 3 TO 1 


4 DESCRIPTIVE NOTES (Type of report and Inclusive dates) 


Scientific. Final. 
B39 AUTHORIS) (First names, middie inttist, last neme) 


Arthur J. Wennerstrom, Richard M. Hearsey 


7a. TOTAL NO. OF PAGES 7b, NO. OF REFS 1 
177 9 : 


Sa, ORIGINATOR’S REPORT NUMBERIS) 


98. OTHER REPORT NOt(S! (Any other numbers that may be assigned i 
thia report) 
ARL 71-0061 


Approved for public release; distribution unlimited. 


6- REPORT OATE 


March 1971 


8a. CONTRACT OR GRANT HO. 
In-house research 
B oS. PROJECT NO. 


B 7065-09 
* DoD Eiement 61102F 


|. DoD Subelement 681307 


R10 DISTRIBUTION STATEMENT 


Pit. SUPPLEMENTARY NOTES 


TECH OTHER 


12. SPONSORING MILITARY ACTIVITY 


Fluid Dynamics Facilities Rsch Lab 
ARL/LF, Wright-Patterson AFB, OH 


This report describes in detail the aerodynamic design of a super- 

fs sonic axial compressor stage. The principal design-point characteris- 

tics of the stage are a corrected tip speed of 1600 ft/sec, an inlet 

, hub/tip radius ratio of 0.75, a total pressure ratio of 3.0, and an 

isentropic efficiency of 82%. Four features distinguish this stage 

from other reported stages. A new type of rotor airfoil is employed. 

| The stator leading edges are swept back from both walls toward mid- 

| passage. Unusually large and variable fillet radii blend blades with 
piatforms. Also, a new and precise technique was used to determine 

{ Cartesian manufacturing coordinates for the airfoils, aerodynamically 

defined on streamsurfaces. The preliminary design employed a technique 

resulting in equilibrium radial distributions of loss coefficient and 

1 flow angle which are fully consistent with relative Mach numbers and 

diffusion factors for each blade row and on each streamsurface accord- 

ing to a prescribed loss model. The detail design was accomplished 

using computing stations internal as well as external to both blade 

rows and attempted to optimize the axial distribution of static 

f pressure. 


DD "1473, ae UNCLASSIFIED 


Security Classification 


¢ 4” 


satan I 
4 my 


Cae 


5 


mone 7 me 
eae 
a BASES hab Sifted 


aa a AES oust 


ae 


Security Classification 


REY wWORO!S 
Proce] wt [rove] wz [rove] wr | 


AXIAL COMPRESSOR ; 
SUPERSONIC COMPRESSOR 

GAS TURBINE 

TURBINE ENGINE 


UNCLASSIFIED 


Security Classification 


Se eeeemeaitl ainesentl 


err 5 
a eR aE 
pion fi. SO ee 


ARL 71-0061 


THE DESIGN OF AN AXIAL COMPRESSOR STAGE 
FOR A TOTAL PRESSURE RATIO OF 3 TO 1 


A. J, WENNERSTROM 
R. M. HEARSEY 


FLUID DYNAMICS FACILITIES RESEARCH LABORATORY 


MARCH 1971 


Approved for public release; distribution unlimited. 


AEROSPACE RESEARCH LABORATORIES 
AIR FORCE SYSTEMS COMMAND 
UNITED STATES AIR FORCE 
WRIGHT-PATTERSON AIR FORCE BASE, OHIO 


fee e™ sce dite 
2 


oes ve 


eaten te 


= Bin techiatt ts a eS Peta geen let 


vbr 


en <3 re ee 
FR Th, dg oy ae amet he: 
Fog ent BOP LTRs 


BS 


Bf 


a % 
et 


% Scat 
peat Hote Fee: 


& 


TER, 
py ae os Se 
Lae, wh ecb BBA 


. 
etre 
pels ae 


at 
BEE 


3 


a One 


S MEMOS oar 
NESS earn ee 


ee 
E + er ant = pee? 
Ecce tbe Re WIP SS ERLE 


<€ 


fs TE 
we et ates 


FOREWORD 


This report was prepared by Mr. Kichurd M. nearsey, Visiting 
Research Associate of the Ohio State University Research Foundation, 
and Dr. Arthur J. Wennerstrom of the Fluid Dynamics Facilities Re- 


search Laboratory, Aerospace Research Laboratories, Wright-Patterson 
Air Force Base, Ohio. 


The report presents results from a portion of the effort of the 
Fluid Machinery Research Group, supervised by Dr. Arthur J. Wennerstrom 
and was conducted under Work Unit 09 of Project 7065, "Aerospace Simu- 


lation Techniques Research," under the over-all direction of Mr, bhimer 
G. Johnson. 
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ABSTRACT 


This report describes in detail the aerodynamic design of a 
supersonic axial compressor stage. The principal design-point char- 
acteristics of the stage are a corrected tip speed of 1600 ft/sec,. 
an inlet hub/tip radius ratio of 0.75, a total pressure ratio of 3.0, 
and an isentropic efficiency of 82%. Four features -Aistinguish this 
stage from other reported stages. A new type of rotor airfoil is em- 
ployed. The stator leading edges are swept back from both walis. 
toward mid-passage. Unusually large and variable fillet radii blend 
blades with platforms. f&lso, a new and precise technique was used 
to determine Cartesian manufacturing coordinates for the airfoils, 
aerodynamically defined on streamsurfaces. The preliminary design 
employed a technique resulting in equilibrium radial distributions of 
loss coefficient and flow angle which are fully consistent with. rela- 


tive Mach numbers and diffusion factors for each blade row and on each 
streamsurface according to a. prescribed loss model. The detail design 
was accomplished using computing stations. internal as well as. external 
to both blade rows and attempted to optimize the axial distribution of 


static pressure. 
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SECTION I 
INTRODUCTION 


A major goal associated with research into and development of 
aircraft~related turbomachinery is reduction of the size, weight and 
cost of the compressor component. Reductions of size and weight 
would benefit many types of equipment, but most particularly aircraft 
turbine engines, gas-turbine powered ground support equipment and 
vehicles, and other air-transportable turbocompressor systems. Such 
reductions may be achieved through increases in specific performance 
levels such as pressure ratio per stage and flow per unit frontal 
area. Cost reductions accrue through two means. First, a reduced 
size and/or number of fabricated parts leads to lower component costs 
to meet a particular performance objective. Second, lower component 
weight allows higher payloads for associated vehicle systems. In- 
creases in specific performance are achieved through increases in 
rotational speed and aerodynamic loading levels, normaliy at the ex- 
pense of thermodynamic efficiency. This report addresses itself to 
the problem of raising thermodynamic efficiency in axial compressor 
Stages designed for high Mach numbers and high pressure ratios. 


The design of a heavily-loaded, supersonic axial compressor 
stage is described. The stage target performance calis for an cver- 
all total pressure radio of 3:1, at an isentropic efficiency (based 
upon stagnation conditions) of 6.82. The rotor corrected tip speed 
is 1600 feet per second, with an inlet hub-to-tip radius ratio of 
0.75:1. This performance level is beyond the current state-of-the- 
art. Accordingly, the intent is to produce a stage having the best 
obtainable performance, and hence with the maximum potential for 
development to the objective performance level. 


Axial turbomachine design is generally considered to consist of 
first, the execution of an aerodynamic design procedure and, second, 
the determination of the geometry of blading that is compatible with 
the aerodynamic design and with pertinent mechanical restraints. 


These two steps were followed in an integrated manner for this design. 


Initially, an aerodynamic design was produced which satisfied the 
design »bjectives and criteria, calculations being made at the blade 
leading and trailing edges. This calculation included the use of 
selected loss models for the rotor and stator blades. Hence the 


results included relative inlet and outlet flow angles and total pres- 


sure loss coefficients for the rotor and stator that were consistent 
with the design assumptions. Subsequently, further aerodynamic cal- 


culations were made at locations both exterior to and within the blade 


rows to evaluate possible blading geometries and details of annulus 
wall contours within the blading. 


The aerodynamic calculations were made using an axisymmetric 
flow analysis computed by the streamline curvature technique. This 
Method is discussed in Ref 1. (The mechanical evaluation of the de- 
sign was made separately from the aerodynamic design.) 


wits a view fo achieving the previously stated goals, several 
nevei features were ineorporated into the design. A new rotor blade 
seagtilon was employed, a key charaeteristic of which is tnat it con- 
tuins mo licvontinuities in surface curvature. This is intended to 
Minimise toundary layer separation due to abrupt velocity gradient 
changes in the flow at the blade surfaces. Having designed the blade 
sections om streamsurfaces-, an accurate procedure for the determina~- 
tion af Cartesian coordinates for the blade was employed to provide 
data convenicat ta the manufacturing process. A full description of 
the Liade section and the method of determining manufacturing coordi-~ 
nate data is given in Ref 2. A swept leading edge configuration was 
included in the stator blade design. It is hoped that this will re- 
duce the sensitivity of the stator blades to incidence at elevated 
Mach numbers. 


Subsequent sections of the report describe the various aspects 
and phases of the design. The criteria incorporated into the design 
are presented followed by details of the aerodynamic evaluation 
method, The results of the various design calculations are described, 
and then the final stage configuration is given. Finally, the antici- 
pated performance of the compressor is discussed, 


SECTION If 
DESIGN CRITERIA 


1. LOADING AND LOSS DISTRIBUTIONS 
a. Qver-all Averaged Quantities 


The target performance for the stage (outlined previously) 
together with some leading parameters specified at the outset of the 


design process, are as follows: 


Total pressure ratio 3.0:1 


Isentropic efficiency (total-to-total) 0.82 


Inlet hub/tip radius ratio 0.75 

Rotor tip speed (corrected to SLS) 1600 ft/sec 
Inlet Mach number 0.55 

Exit Mach number 0.5 

Exit swirl zero 


An 18-inch diameter was selected for the compressor, resulting in a 
design flow-rate (corrected to SLS) of 30.0 lbs/sec. (This dimension 
resulted from consideration of power requirements for testing, as an 
upper bound, and a practical minimum size for instrumentation as a 
lower bound. Aerodynamically, only the blockage allowance attributed 
to annulus wall boundary layers is significantly affected.) 


b. Radial Distributions 


Because maximum design-point efficiency is a principal objec- 
tive of this stage, blade element minimum loss values have been used 
for each element, consistent with local values of diffusion and Mach 
number. Losses, expressed as relative total pressure loss coeffi- 
cients, are assumed to be composed of two additive components; a loss 
due to diffusion and a loss due to shock waves (where Mach number is 
high enough). Then, the over-all blade element loss coefficients for 


rotor and stator are defined. 
® = oO, + O 
(1) 
is the over-all relative total pressure loss coefficient 


is the contribution due to diffusion losses 
is the contribution due to shock losses. 
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Yhe diffucion less component model is based upon the well-known 
WACA Correlation presented in Ref 3. A total pressure loss param- 
eler and a diffusion factor are defined by 


Oa cos Son 


P= a ean (2a) 
and : 
foe ee ee eee (2b) 
ae 2C, po 
where P is the total pressure loss parameter 
os is the relative outiet flow angle 
o is the cascade solidity 
D is the diffusion factor 
c. is the relative outlet velocity 
Cos is the relative inlet velocity 
pan is the relative inlet whirl velocity 
Cen is the relative outlet whirl velocity. 


The relationship that was assumed between the total pressure loss 
parameter and the diffusion factor is shown in Fig 1, This is an 
extrapolation to higher diffusion factors of the data given in Fig 203 
of Ref 3. 


The shock loss component is determined from the loss across a 
normal shock, the strength of which is a function of the blade ele- 
ment relative inlet Mach number. Different blade profiles were 
selected for the rotor and stator (as discussed in the following 
section) and hence two different relationships for the shock Mach 
number as a function of the inlet Mach number were used. 


the profile selected for the rotor blade has almost no cam- 
ber upstream of the point where the main passage shock will probably 
impinge upon the blade suction surface so that, from the viewpoint 
of a two-dimensional cascade, the Mach number should be unchanged 
relative to the inlet value. Compounded with this are effects in the 
meridional (axial-radial) plane. ‘The annulus walls were assumed (at 
this phase of the design) to be converged in the region of the rotor 
leading edge co that significant contraction and hence diffusion of 
the supersonic relative flow would occur. The average Mach number 
immediately upstream of the shock is assumed to be given by 


ANRITSU ILL Nee AE MERE a rt Re 


An ee, 


OVI Sa re 


POO A MR tn Al ERE IO Ea a BIER URE TPEN BO te 
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M, = 1.0 + 0.6667 (Me - 1.0) (3) 
where M. is the Mach number immediately upstream of the shock 
Me is the blade section relative inlet Mach number. 


Then the relative total pressure loss coefficient component due to 
shock losses is given by 


y+ lyme af ye } 2 
Ve : ae ee: 7 Y¥-2 ees y-1 
J Gy-2)M2 + 2 2yM2 - (y-2) 
a ae neal (4) 
¥ 
y-1 Qo. abit “Se 
1 - ; en | y-2 


The double circular are profile used for the stator had convex 
curvature throughout the suction surface, and hence, from the view- 
point of a two-dimensional cascade, can be expected to expand a super- 
sonic inlet flow to a yet higher Mach number upstream of the main 
passage shock. Miller, Hartman and Lewis examined shock losses in 
double circular arc rotors and presented (in Ref 4) a simple flow 
model that correlated well with experimental data. A normal shock is 
presumed to lie along a straight line from the leading edge of one 
blade to a shock-impingement point on the adjacent blade suction sur- 
face. The shock-impingement point is fixed by the assumption that 
the shock intersects normally a blade mean~-line drawn in mid-passage. 
The flow Mach number at the shock-impingement point is determined by 
a Prandtl-Meyer expansion from the value at blade inlet through the 
angle subtended by the suction surface upstream of the shock. by 
assuming the shock strength to be given by the mean of the Macn num- 
ber at each end of the shock, the shock loss may be determined. For 
the current design, it was assumed that the result of the two- 
dimensional analysis is further modified by annulus wall effects, 
as in the case of the rotor. The Mach number produced by the above 
analysis is presumed to be modified according to bq (3), the elevated 
mean two-dimensional shock Mach number replacing the blade section 
relative inlet Mach number. Equation (4) is then applied to deter- 
mine the relative total pressure loss coefficient component due to 
shocks in the stator., 


ec. Axial Distributions Within Blade Rows 


The averaged quantities and radial distributions discussed 
above were associated with conditions spanning an entire blade row. 
fhe detailed aerodynamic analyses (described in Section I1IiI) include 
computations made at a number of axial stations internal to both 


St 


wee. 


eee 


blade rows. For this purpose it is necessary to specify the axial 
distribution of total pressure losses along each streamsurface and 
within. each blade row. The quantity chosen to define a distribution 
was the ratio of actual~ to-ideal relative total pressure on a stream- 
surface, referred to conditions at the leading édge of the respective 
blade row. Then the (absolute) total pressure at any station may be 
determined from 


P =P (2) YI : Par 
“n 1 T. ‘ 
: Par (ideal) (5) 
where Pi is the désired total pressure 
i. is the total pressure at blade row inlet 
Ty 
r is thé ratio of total temperatures between 
i the point of interest and blade inlet 
Por 


= is the ratio of actual to ideal relative 
Par(ideal) total pressures at the point of interest 


The value of Par/Pnr( ideal): at the blade outlet was obtained from 
Y 
ae 
Por (ideal) | an stens)) 


mn) 
(6) 
p 


where ye 


P s as 0 (ee 
2r (ideal) = 2 1] A ‘ ‘Y-1 (7) i 


2 r i 
Pir 2 Goin 2 ; : 


where P is the relative total pressure at blade inlet 


is the relative total pressure at the blade outlet 


eer hf are the streamsurface radii at the corresponding 
two points 


is the wheel speed at the blade outlet 


is the relative Mach number at the blade inlet 


ao, is the stagnation speed of sound relative to the 
Pe  plade inlet. 


Equations (5), (6) and (7) are equivalent to-Eqs A2 3 
AG 
respectively, of Ref 3, page 253 gq , and A3, 


Although the losses due to the Passage shock occur abruptly, the 
shock is not perpendicular to the meridional plane of the comnressor 
and hence appears in mean effect as a continuous phenomenon. aArbi- 
trarily, the variation of Por/Pnr(ideal) was made linear between valies 
of unity at the blade leading edge and the determined value for the 
blade trailing edge. For the rotor it was convenient to make it 
linear with axial distance, while for the stator streameurface lense: 
Was utilized, 


2. BLADE PROFILE AND SOLIDITY SELECTION 
a. Rotor Blade 


The rotor blade is required to operate (at design point) at 
Mach numbers varying from 1.22 at the hub to 1.58 at the casing. Ar 
these Mach number levels a prime consideration in the selection of a 
blade profile is the minimization of losses due to shock waves, Thix 
rotor design is based upon the concept of achieving essentially 
shock~-free supersonic diffusion in the forward portion of the rotor 
passage, this terminating with an inevitable strong shock. Then sub- 
sonic diffusion to the exit condition occurs. . . 


Supersonic diffusion implies an area decrease. In this design, 
the area reduction is achieved by a decrease in the compressor annu- 
lus height. The flow angle relative to the rotor decreases by an 
amount sufficiently small that the associated flow area increase in 
the cascade plane is relatively small. It is assumed that the super- . 
sonic diffusion will be achieved by compression waves propagated fror: 
the suction surface of the blade upstream of the shock. The compres- 
sion waves should occur as the flow is deflected away from the suc- 
tion surface, toward which it will tend to move in the presence of | 
the flow area reduction. (A similar hypothesis leads naturall) to 
the use of reverse or negative camber in the forward region of a 
supersonic blade section. A potential advantage of the deslgn metho: 
pursued for this stage is that the amount of positive camber requires 
to achieve a specific outlet angle will cf course be less when it is 
not preceded by reverse camber. A comparison of the use of "annulus 
controlled" area reduction arid the 'S' profile for supersonic rotor 
duty is an area for future research.) 


In the past, the types of profile associated with this design 
philosophy have included multiple-circular-arc profiles, and pro- 
files consisting (on the camber line) of a straight line followed by 
a circular are. For this design a new profile is emploved. This i 
comprised of a polynomial (quartic) camber line, with a thickness 
distribution applied about it that consists of two polynomial (cubic, 
curves, one ahead and one after the point of maximum thickness. Sy 
setting the second derivative of the camber line equal to zere at 
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the leading edge, and to one half of its maximum value at the trailing 
edge, a camber line is produced that is. initially straight and has a 
maximum curvature forward: of the trailing edge. The second deriva- 
tive of the thickness distribution equation is also set equal to zero 
at the leading edge, thus maintaining a small mintegit edge wedge angle 
(assuming a conventional thickness to chord ratio), while preventing 
inverse curvature in tne thickness distribution. At the point of max= 
imum thickness, the section thickness and first, and second derivatives 
of thickness are Set equal for the two thickness equations. The re- 
sult is a blade section with continuous surface curvature throughout. 
Further details of the profile are given in Ref 2. 


It can reasonably be assumed that the terminal shock in the 
rotor passage will be nearly. perpendicular to the relative flow and 
will extend from the leading edge of one blade to the suction surface 
of the adjacent blade. (Correlations of such a flow model with ex- 
perimental data were made :by Miller,, Lewis ‘and Hartmann in Ref 4 as 
mentioned before. ) The design philosophy: outlined previously calls 
for little camber in the region of the blade between. the leading edge 
and shock impingement point. This places a lower limit on solidity 
for the blade profile selected. The leading portion of the blade is 
essentially straight, whatever camber is specified being concentrated 
toward the rear of the blade; As solidity is decreased, the pre~ 
dicted point of shock impingement moves rearward along the blade suc~ — 
tion surface, so that there is a solidity below which the require- 
ment of little camber (in the forward region). is not met. The min- 
imum acceptable solidity indicated by this means is about 1.7 at all 
radii. ‘ 


The diffusion factor loss correlation described above also pre- 
sents a. method of determining solidity. For a given air-angle de- 
sign, the variation of minimum low-speed loss coefficient with solid- 
ity may be obtained. Generally, a shailow minimum exists at a 
solidity which increases with diffusion factor. Using this method, 
"optimum" solidities of 1.3, 1.4 and 1.0 at hub, mid and tip respect- 
ively may be calculated. (These figures are consistent with the 
velocity triangle data produced by the final run in the first phase 
of the design. Similar figures were calculated from preliminary re-~ 
sults.) Intuitively, these figures appear to be somewhat low, and 
may represent a minimum feasible solidity for low~speed flow. 


Mechanical considerations also affect the choice of rotor solid- 
ity distribution. The requirement. of a high solidity at the tip sec- 
tion conflicts with the requirement to keep the ratio of the cross- 
sectional areas of the casing and-hub sections down to such a value 
that unacceptable centrifugal stresses are not generated in the rotor 
blade. 
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Taking into. account the three factors described above, the fol- 
lowing distribution of rotor blade solidity was derived. 


Section Solidity 
Hub 2.272 
Mean 1.937 
Casing 1.861 


Evidence presented in Ref 5 indicates that no loss attributable 
to trailing-edge thickness was observed for compressor cascades 
having trailing-edge thicknesses up to about one-third of the maxi- 


mum blade-element thickness. Therefore, for this design, rotor trail- 


ing edges were simply truncated at one-third of maximum blade-élement 
thickness in-order to compensate to some extent for the annulus con- 
vergence required between rotor and stator in order to achieve smooth 
wall contours. 


b. Stator Blade 


The stator blade operates (at désign point) at relative inlet 
Mach numbers between 1.0 and 1.1. The double circular arc profile 
has performed well at this Mach number level. It constitutes a good 
compromise between expanding the incoming flow to a yet higher Mach 
humber, and maintaining a sufficiently large throat width to pass the 
design flow without choking. It was therefore selected for the stator 
design. 


The same considerations applied to the rotor to derive solid- 
ities may again be used to obtain solidities for the stator. From. a 
two-dimensional viewpoint, emphasis should be placed on minimizing 
the expansion of the incoming supersonic flow upstream of the ter- 
minal shock. As the blade suction surface is a circular are, the ex- 
pansion is continuously decreased as the solidity is increased. How- 
ever, the modest inlet Mach number precludes the assumption that an 
extreme solidity is desirable. The diffusion factor loss coefficient 
correlation may again be used to determine optimum solidities, and 
indicates values varying from 1.2 at the hub to 1.9 at the casing. 


The final stator solidity distribution was influenced consider- 
ably by the use of the swept leading edge, described later. However, 
the concepts described above guided the selection of the general 
solidity level. Final stator solidities are as follows> 


Section Solidity 
Hub Pa3oS 
Mean 1.803 
Casing 2.145 
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3. DEVIATION ANGLES 
a. Rotor Blade 


In orcer to relate blade angles to design air angles it is 
necessary to have a knowledge of the flow deviation angles that. will 
occur. This presents a problem for the rotor design as the profile 
type selected has not been previously investigated and hence a cor- 
relation of deviation with camber, outlet angle, solidity, and so 
fourth, is not availablé. Fortunately, at the loading level for 
which the rotor is designed, a Small error in deviation angle is not 
crucial. An important result from the testing of the stage will be q 
the détermination of the actual rotor deviation angles achieved. ; 


A form of Carter's rule was used. to determine rotor blade 
deviation angles. Deviation is related to cascade geometry by 


gd. ge Aare (8) 
where 6 is: the deviation angle 
6 is the blade section camber angle 
0 is the cascade solidity 
m is a function of the blade. section stagger angle. 


Figure 2 shows the relationship that was assumed between m and the ; 
stagger angle. (Also shown are Carter's curves for conventional 7 
tbliades having their points of maximum camber at 0.4 and 0.5 of the 
chord. These are taken from Fig 160, Ref 3.) 


In order to perform detail design calculations on. proposed rotor E 
configurations, it is necessary to have a means of determining mean q 
relative flow angles within the rotor blade, This was accomplished 
by. uSing a generalized relationship for the deviation angle (from 
the local blade camber line direction) as a function of distance 
along the section meridionally=proj.ected chord and the final devia- 
tion angle. The relationship assumed is shown in Fig 3. The rapid 
increase of deviation angle near the blade exit arises from consid- _ 
eration of the physical requirements necessary to satisfy the "Kutta" 
condition at the trailing edge. 


b. Stator Blade 
Tne same requirements of Knowing final and intermediate devia- 


tion exist for the stator blade as for the rotor blade. However, in 
this case the double circular are profile was specified. 
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Final deviation angles were calctiated from Carter's rule, 
which was described above in connection with the rotor blade. In 
this case the m/stagger curve for conventional blades having maxi- 
mum camber at mid-chord was used. ‘Carter's rule has been widely and. 
successfully used, for double circular arc blade sections. 


Within the stator blade, the same generalized reiationship for 
local deviation angle was assumed as was described above for the 
rotor blade. 


4. ASPECT RATIOS, FILLET RADII, 
AND STATOR LEADING EDGE FORM 


Some novel features of the design are related to attempts to 
minimize the detrimental effects of secondary flows within the com- 
‘pressor. Conclusions. drawn from data presented in Ref 6 inspire 
these attempts. A supersonic radial compressor diffuser of unusually 
high. pérformance is described. The diffuser is comprised of a series 
of uniformly distributed eircular cylindric passages in a radial 
plane, all tangent to a circle in that plane. They mutually inter- 
Sect at the inner radius: :or inlet of the diffuser, forming a series 
‘of sharp, elliptical leading edges. The significant features which 
may have contributed to the performance of the diffuser are believed 
to be the use of a circular (corner-free) passage, and the "swept" 
leading édge configuration, 


The circular passage concept was applied to both rotor and 
stator. It was incorporated by. choosing the number of blades so that 
the flow passages are approximately square when viewed normal to the 
flow near the blade row exits. Also large fillet radii were applied 
at the rotor hub, and the stator hub and casing. 


The swept leading edge concept has been applied to the stator 
blade. Swept wing theory indicates that section performance is re- 
lated to the normal component of the incident Mach number. High 
Mach: number operation of compressor: blade séctions is characterized 
‘by: relatively large loss pénalties for small variations in incidence 
angle away from the optimum, The hub and the casing are regions 
where the stator design incidence is most likely to be violated. 
Accordingly, the stator leading edge at the hub and casing has been 
swept forward to such an angle that the normal Mach number is approx- 
imately 0.4, A simple parabolic form connects the two extremities, 
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SECTION If! 
AERODYNAMIC CALCULATION METHOD 


1. AXISYMMETRIC FLOW ANALYSTS 


The principal means used to incorporate the design criteria into 
the stage design was an axisymmetric flow analysis. The most impor- 
tant assumpLions made are that the flow is axisymmetric, and that 
there is no transport of mass or energy across streamsurfaces in the 
flow. The fluid (atr) is assumed to be a perfect gas. Briefly, the 
calculation consists of the following elements: 


(1) A number of computing stations are located in the flow 


(2) The lhecations of a number of axisymmetric streamsurfaces 
are estimated at each computing station 


(3) The continuity, momentum, and energy equations are simul- 
taneously solved (iteratively) at each station in turn. The conti- 
nuity eauation is satisfied in an integrated sense, that is, the 
svecified flow-rate at each station is maintained. The energy equa- 
tion is satisfied by one of essentially three means. At a station 
following a blade-free space, the enthalpy, entropy and angular 
momentum are constant along streamsurfaces from the preceding station. 
At a station within or immediately following a blade row, there may 
be specified (for each streamsurface) either the work input down- 
stream of a preceding station, or the flow angle relative to the 
blade. In both cases. the angular momentum and enthalpy change are 
established, directly or indirectly. A number of means of specifying 
the entropy rise on each streamsurface exist. For the current purpose, 
two alternative specifications were sufficient; an inlet dynamic head 
total pressure loss coefficient, or the ratio of actual to ideal 
relative total pressure. The solution to the momentum equation, 
which may be considered to be the principal equation, yields the var- 
fation of velocity, and hence all other undetermined parameters of 
the flow, along the computing station. 


(4) The estimated streamsurface pattern used to obtain the 
solution described in (3) is refined to more nearly satisfy continuity 
on a detailed basis. That is, using the mass flux distributions de- 
rived in (3), the streamsurface location estimates are revised to 
maintain a constant proportion of the total flow in each streamtube. 


(5) Tne procedure is re-entered at (3) to obtain an improved 
SOlution to the system of equations. The new solution differs from 
that previously determined because it is a function of the assumed 
streamsurface pattern. This is repeated until the desired accuracy 
is achieved. 
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where 


The momentum equation used for the above calculations is 


cos(y+¢) tana d(r tan a) sin(y+o) dC. 


= 2 cos*aC* |—__——. - —_———— + ————- fae + 
a a is de C. m 
(9) 
2. | dH dS 
as ee 
2 cas aE 4 


C. ais the meridional velocity 

sa is the direction of the computing station 

a is the whirl angle, defined tan a = Co/Cin 

C, is the tangential velocity 

Y is the angle made by the computing station with the radial 
direction, positive values indicating an. increase in 
radius with axial distance 

0) is the streamsurface slope angle 

r, ais the radius of curvature of the streamsurface 

m is the meridional streamline direction 


H is the enthalpy 


S is. the entropy 


et 


is the static temperature. 


The continuity equation is 


case 
are C,, cos @ w dA (10) 


hub 


where W is the flow-rate 


W is the specific weight 


A is flow area normal to the axis. 
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The energy equation is either 


y = constant along stream- 
lines between blade rows. (11). 

specified value within 

T = or immediately following (12) 

a rotor row 
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when relative flow angle is specified. In Eq (13), fe 
T is total temperature | 
U is blade speed ; 
ih 
g is the acceleration due to gravity 
J is Joule's equivalent 4 
: 4 
C is specific heat at constant pressure 
n indicates the location where the temperature is derived 
indicates ‘blade inlet 
r indicates relative conditions. 
Equations (5), (6) and (7) (given previously) relate total pressure 
to total temperature. 
2. IYERATION PROCEDURE FOR CONSISTENT LOSSES 
One of the problems faced by the designer of a compressor (or i 
turbine) is how to specify losses in such a manner that they are com- ; 
pletely consistent with local values of Mach number and diffusion 
along every streamsurface. When using a streamline-curvature calcu- 
lation technique, it is convenient to solve this problem iteratively. 
The general procedure followed was: 


(1) A complete aerodynamic solution was obtained using initially 
estimated loss distributions (which could be zero loss throughout) 


(2) Mach numbers and diffusion factors were calculated on each 


streamsurface for each blade row and new loss coefficients were ob- 
tained using Eqs (1), (2), (3) and (4) ss 
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(3) A new aerodynamic solution was obtained using the loss co- 
efficients obtained in step (2). 


Steps (2) and (3) were repeated until the change in loss coefficients 
fell within a prescribed tolerance. For these calculations,. comput- 
ing stations within blade rows were not employed; only those corres- 
ponding to blade-row leading and trailing edge planes and elsewhere 
in the compressor were used. The radial distribution of any one of 
three parameters could be defined in order to control stage perform- 
ance. The three options were to specify: 


(1) Stage total pressure ratio at the stator trailing edge plane, 
(2) Total temperature rise at the rotor trailing edge plane, or 


(3) Relative flow angle at the rotor trailing edge plane. 


3. ALTERNATIVE FORMULATION 
OF MOMENTUM EQUATION 


Some difficulty was experienced in selecting a radial distribu- 
tion of work (total temperature or pressure rise) that would yield a 
satisfactory velocity profile at the rotor exit when making the com- 
putations described in the previous section. An alternative formu- 
lation of the momentum equation was therefore derived in which the 
total temperature rise is the dependent variable, and the velocity 
profile is specified. The equation to be solved is formed by combin- 
ing Eqs (9), (13) and (5), and using the fact that for a perfect gas 
we may write 


ea: eee t ( at 1 a 
de” ae GC a > 7 (oe, oe Pn de (14) 


Algebraic manipulation yields the following result 
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Thus the total temperature gradient is a function of the meridional 
velocity gradient, the gradient of actual to ideal relative total 
pressure (a function of the losses), and other quantities. By choosing 
@ value of total temperature at one radius, the values at all other 
radii may be calculated. The starting point value is adjusted. so that 
together with the iteratively determined losses, the desired mean. 
stage total pressure ratio is achieved. This is further described in 
Ref 7. 


4, ANNULUS WALL BOUNDARY LAYER DETERMINATION 

A constant phenomenon in axial compressors is the build-up of 
boundary layers of significant thickness upon the annulus walls. The 
foklowing simple calculation was included in the computing scheme used 
for the stage design. 

Jansen presents a method of estimating the blockage due to annu- 
lus wall boundary layers in Ref 8. Boundary layer displacement thick- 
nesses are specified at a location upstream of the compressor: proper. 


Here they will generally be negligible. Then the momentum thickness 
ef each boundary layer is obtained from 


Mn 70.8 
6 = 8, + 0.006 C “sal f cg aa 
n mn m™ : (16) 


- m, 


where 6 is momentum thickness. (Note that units of feet and seconds 
are assumed for the dimensional quantities.) 


The shape factor is obtained from 


e : a 
H. = 30 ————— + 1.5 (17) 


subject to the restraints 1.1 < H < 2.2 where H is the shape factor. 


The displacement thickness is given by 


n nn (18): 
Where 6 is displacement thickness. 
The blockage due to the two boundary layers is incorporated into 


the calculation by locating the outermost streamsurfaces away from the 
annulus walls by the amount of the displacement thickness. 
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5. USE OF STATIC PRESSURE DISTRIBUTIONS ‘TO 
OPTIMIZE INCIDENCE ANGLES AND ANNULUS GEOMETRY 


Given the results of the iterative loss reestimation procedure, 
it is possible to determine feasible blade geometries by assuming 
incidence angle variations for rotor and stator. The calculation 
method employed includes the calculation of conditions, on an axisym- 
metric basis, at points within the blade rows. Hence, within the 
limits of the assumptions, it is possible to determine the variation 
of any parameter through. the blade rows. A rational method of evalu- 
ating various designs was required to enable incidence variations and 
annulus configurations to be optimized. 


The parameter selected for prime consideration was the static 
pressure. The basic concept employed in the optimization process was 
that the static pressure should rise in a smooth manner with minimum 
slope. ‘This minimum was limited by the requirement that the rate of 
increase of static pressure with flow-path length should fall smoothly 
to zéro at the blade trailing edge. The validity of this approach is 
debatable for the rotor blade sections, which operate transonically 
with a strong shock in the flow. However, the shock is approximately 
normal to the relative flow direction and hence, when viewed in the 
meridional ‘plane and considered in mean effect, appears not as a dis- 
continuity but as a region. through which conditions change in an 
apparently continuous manner. 


An implicit result of the intra-blade analysis is a check upon 
the maximum or choking flow of the ‘blade row. 


6. COMPUTER PROGRAM 


A computer program to perform the calculations described above 
was created by modifying the program presented in Ref 9, which de~ 
seribes a program for the analysis of non-axisymmetric flows in axial 
compressors, The program is considerably more complicated than the 
axisymmetric analysis outlined above requires, and extensive simplifi- 
cations were made to the deck. The principal changes made were as 


follows. 


(1) The system of equations solved was modified to reflect the 
assumption of axial symmetry. (The equations, in the modified form, 
were presented previously (Eqs (9), (10), (11), (12) and (13) 


(2) The loss estimation procedure described previously was pro- 
grammed and incorporated into an overall iteration procedure. 


(3) The method of solution of the momentum equation in which the 


dependent variable is the total temperature gradient (Eq (15)) was 
programmed as an alternative calculation at the rotor exit. 
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(4) One novel feature incorporated into the program is that : 
the station lean angle (y) may vary along 4 computing station. This: 
provision was included in order to be able to locate calculating 
stations along the leading edge of the stator and within the stator 
blade row at constant fractions of the projected blade chord. 1 4 
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SECTION IV 
RESULTS FROM DESIGN CALCULATIONS 


1. ITERATIVE LOSS REESTIMATION PROCEDURE 
a, Path to a Solution 


As mentioned in the preyious section, several different input 
options to the calculation scheme (i.e., computer program) were avail- 
able that would determine, directly cor indirectly, the rotor work dis- 
tribution. The criteria for selection of an option for the first 
phase of the design are that the mass-averaged tatal pressure ratic of 
the stage should be readily controllable, and that satisfactory radial 
distributions of the significant parameters of Loth rator and stator 
should be easily obtained. The latter criterion is closely allied to 
the rotor exit meridional velocity profile. 


Specification of the distribution of either total temperature 
rise across the rotor or relative outlet flow angle from the rotor 
satisfies neither criteria. Specification of the distribution of 
rotor (or stage) outlet total pressure gives relatively goad control. 
For this stage design, the specification of a nondimensionalized rotor 
outlet meridional velocity profile together with a mass-averaged total 


pressure ratio was selected, this being thought to be the most direct 


method of achieving the desired result. In faet, the total pressure 
distribution is probably as easy to manipulate. 


Having achieved a satisfactory aerodynamic result by the above 
means, it remains to ensure that the blades implied by the forcgoing 
calculations are mechanically acceptable, especially in the case of 
the rotor blade. The rotor relative flow angles must vary along the 
length of the blade in such a manner that no severe blade twicts occur. 
In order to achieve this for the current desigr, aerodynamic analyses 
were made for specified relative rotor outlet angle distributions 
which were produced by "smoothing" the distributions determined in the 
preceding calculations. This resulted in a stage mass-averazed total 
pressure ratio little different from the defined value of 3.9:1, and 
mechanically acceptable blade shapes. 


b. Design-Point Conditions 


An important result of the calculations described above was 
the over-all stage efficiency produced by the loss model discussed 
earlier. Actually, it was originally anticipated that a more op- 
timistic diffusion-loss model would be required to produce the objec- 
tive stage efficiency. However, this was achieved with the loss-model 
described, which is believed to be quite realistic as far as diffusion 
loss is concerned, if somewhat optimistic with respect to shock losses. 
Some degree of optimism may also be involved in the diffusion loss 
model inasmuch as the elevated tip loss incorporated in the NACA model 
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“ued te waxiniging the static pressure rise across the rotor while 
“it qaintnainius satisfactory rotor diffusion levels. This resulted 
Vt TA ty fivwher static pressure rise across the rotor than across 
Tees 'atory but the mean meridional velocity rises over 100 feet per 
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2. INPRA-BLADE FLOW ANALYSIS 


The objective of the intr ~blade flow analysis was to determine 
details of blade and annulus geometries that were consistent with the 
results obtained from the iterative loss re~estimation procedure, and 
would also satisfy the design criteria stated earlier. The approach 
that was followed consisted of analyzing conditions within the com- 
pressor for possible detailed blade and annulus geometries, and then 
modifying the assumed geometries until the desired results were 
achieved. A somewhat lengthy trial-and-error process is implied and 
was indeed necessary. 


The input data to the computer program consisted of the annulus 
geometry, the performance of the blading, the flow-rate, and the rota- 
tional speed of the machine. The annulus geometry is, of course, 
specified by giving the radius of the hub and of the casing at each 
computing station. Hence this was adjusted by changing these radii 
within the blade-rows, the values elsewhere being those established by 
the results of the previous design phase. The flow boundaries used 
were the inner limits of the boundary layer displacement thicknesses, 
previously established. A linear variation in displacement thickness 
with length was assumed within the blade rows. The annulus geometry 
finally derived and the computing stations used for the calculations 
are shown in Fig 14. The performance of the blading is specified by 
giving, at each computing station within or immediately following a 
blade row, the radial variations of the relative flow angle, the ratio 
of actual to ideal relative total pressure, and the blockage due to 
the blades. The relative flow angle was obtained from the local cam- 
verline direction, the assumed deviation angle at the trailing edge, 
and the assumed variation of deviation angle within the blade. These 
items were all discussed in Section II. The ratio of actual to ideal 
relative total pressure was obtained from the results of the iterative 
loss re~estimation procedure and Eq's (6) and (7). The blockage due 
to blades is the ratio of blocked to total circumference. 


The computer output obtained for the analysis of the configura- 
tion that was finally selected is shown on the following pages. Com- 
patibility between the intra-blade calculations and the inter-~blade 
calculations performed with the iterative loss re-estimaticn procedure 
is shown by plots of the meridional velocity profile at the rotor 
leading and trailing edges (Fig 15), and the stator leading and trail- 
ing edges (Fig 16). Results from both calculations are displayed, and 
it is seen that differences are quite small, except at the stator 
leading edge. At this computing station in particular, the streamline 
characteristics are computed to be significantly different in the two 
cases. This is due to details of the annulus wall contours which are 
not apparent to the inter-blade calculations. 


The static pressure distribution through the stage (which was 
considered to be a key criterion) is shown in Fig 17. The established 
design criteria regarding the static pressure rise through the blade 
rows is reasonably well satisfied. The resulting incidence and devia-~ 


Pe a 


tion angle variations for the rovor and stator are shown in Figs 18 
ang 19. 
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SECTION V 
FINAL STAGE CONFIGURATION 


1. DESIGN POINT 


The following figures are consistent with the data given in. 
Section ITI: 


Design speed (corrected to SLS) 20371.4 rpm 
Design flow (corrected to SLS) 30.0 lbs/sec 
Design total pressure ratio 3.0 $2 


2. ANNULUS GEOMETRY 


Figure 20 shows the flowpath ineluding inlet bell mouth, compres~ 
sor proper, and exhaust that will be incorporated into the compressor 
test facility. Details of the inner wall (hub) of the compressor 
annulus are shown in Fig 21; Fig 22 shows details of the compressor 
outer wall (casing). In each of thése figures the origin of the axial 
coordinates is the same as was used for the aerodynamic design calcu- 
lations; that is, the rotor leading edge. 


3. ROTOR GEOMETRY 
a. Number of Blades 


The rotor contains 30 blades. 


b. Blade Form 


The rotor blade design was produced by determining profiles on 
15 streamsurfaces, and then stacking the sections to complete the 
blade definition. In order to create data convenient for the ‘manufac- 
turing process, coordinates of plane sections through the blade per- 
pendicular to the stack axis were interpolated (and extrapolated) by 
a "spline-curve" method. A brief description of the (streamsurface) 
airfoil sections was given in Section II of this report; Ref 2 gives 
a full description of the section, the stacking procedure, and the 
method of interpolating (or extrapolating) the "manufacturing" section 
data. A computer program to perform the necessary calculations was 
also presented in Ref 2. For this design, it was convenient to use a 
slightly modified form of this program. The program was amended so 
that the entire blade construction was performed a specified number of 
times. Input data to each pass (after the first) was based on the 
results of the preceding pass so that the axially-projected chord of 
each section would equal the desired value of 2.0 (inches), and so that 
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the blade leading edge, at 211 radii, would lie at the desired loca-~ 
tion, that is x = 0. Changes in the original input data were made to 
the stacking point location, the meridionally~-projected secticn chord 
lengths, and the x-direction stacking offsets of the sections. Thus 
the resulting blade exactly matched the allocated space in the aero- 
dynamic calculations. These iterations could have been made "ty hand" 
using the program as published in Ref 2; the modified program enabled 
the desired result to be obtained from one run. Five iterations were 
made;. in. fact an acceptable result for ali practical purposes was 
achieved after two or three iterations. 


Shown on following pages is computer program output for the rotor 
blade design. All dimensions are in inches. First appear sundry con- 
stants and a definition of the 15 streamsurfaces. The streamsurfaces 
are defined at eight axial locations which coincide with eight of the 
computing stations used for the aerodynamic design caiculations. The 
origin for the axial locations of the stations is the same as was used 
for the aerodynamic analyses. The input data printout is completed 
with a table defining the geometry of each section. A detailed de- 
scription of the significance of each input data item is given in 
Ref 2, but it should be noted that, in accord with the program modifi- 
cation described above, the stacking point location, the section 
meridionally-projected chord lengths, and the x-direction offsets are 
first-estimates only. Next are shown details of the 15 streamsurface 
sections. Only the "normalized" data has been reproduced; the equiva- 
lent dimensional data would be derived by scaling the non-dimensional 
quantities by the meridional chord of the section (or the appropriate 
power thereof). Finally, details are shown of 11 manufacturing sec- 
tions through the blade. These plane sections perpendicular to the 
stack axis are spaced 1/4-inch apart, and extend slightly beyond the 
blade in both directions. The "Z" coordinate is measured along the 
stack axis from the machine axis. The origin for the section coordi- 
nates is the stack axis. The "X" direction is parallel to the 
machine axis, and the coordinate increases in the direction of flow. 
The "Y" direction is perpendicular to the "X" direction, and the "Y" 
coordinate decreases in the direction of rotation. "XS" and "YS" de- 
fine the suction surface of the section, and "XP" and "YP" define the 
pressure surface, "XSEMI" and "YSEMI" define the leading edge radius. 
The trailing edge is a straight line joining the pressure and sucticn 
surfaces. The data reproduced shows 50 points per blade surface; for 
manufacturing purposes the program was run with 120 points per sur- 
face specified. 


Figure 23 shows superimposed piots of developed streamsurface 
sections 1, 3, 5... 15. Figure 24 shows similar views of manu- 
facturing sections 1, 3, 5... .d11. Extrapolation of the blade to 
planes beyond the hub and casing causes the larger changes in section 
along the blade that are seen in Fig 24, relative to Fig 23. 
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| ec. Location of Stack Axis i 
The rotor stack axis is located at an axial coordinate of 4 
0.9791 inches, measured from the same origin as was used to define : 

the annulus geometry. 


ad. Root Fillet 


Between points 3/4-inch in from the leading and trailing edges ‘ 
of the blade, the root fillet is a 1/4-inch radius, on both sides of : 
the blade. The fillet is smoothly decreased to 1/16-inch radius at ‘ 
the blade edges. ; 
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Rh, STATOR GEOMETRY 
a. Number of Blades 


The stator contains 49 blades. 


b. Blade Form 


The stator blade design was produced by the same process as was : 
used for the rotor design, but with two differences. First (as de- a 
seribed in Section II of this report), the double circular arc profile Ft 
was selected for the streamsurface sections. Second, the sections 
were stacked at the trailing edge to complete the blade definition. 


The 15 streamsurfaces on which the double circular are profiles 
were designed are defined by the results of the aerodynamic analysis 
given in Section IV.2 of this report. The stator occupies stations 
11 through 15. Details of each section are given in the following 
table. For each section, the ratio of maximum thickness to chord is 

-O4, and the edge radii are .005 inch. 


Camber Stagger 
section Angle Angle Chord 
No. (degrees) (degrees ) (inches) 
1 (Hub) 54.39 18.790 2.349 
2 54.59 18.54 2.240 
3 54.77 18.38 2.143 
4 55.03 18.26 2.060 
5 55. 33 18.16 1.990 
6 55.59 18.08 1.934 
7 55.85 18.03 1.895 
8 56.13 18.02 1.873 
9 56.39 18.05 1.869 
10 56.68 18.13 1.885 
11 56.99 18.28 1.923 
12 57.49 18.54 1.984 
13 58.35 18.96 2.072 
14 60.38 19.31 2.197 
15 (Casing) 62.35 20.73 2.365 
123 
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The generation of the blade geometry including interpolation 
(and extrapolation) of "manufacturing plane" data was accomplished 
using a computer program similar to that used for the rotor blade, 
and output from the program is shown on following pages. The coordi- 
nates of 10 "manufacturing" plane sections through the blade perpen- 
dicular to the stack axis are given. The sections are spaced 1/8 inch 
apart, and extend slightly beyond the actual blade in both directions. 
The coordinate definitions are the same as given for the rotor blade 
above. 'XP' and 'YP' define the pressure surface of the section, and 
"'XS' and 'YS' define the suction surface. ‘'XSEM1' and 'YSEM1! define 
the trailing edge radius, and 'XSEM2' and 'YSEM2' define the leading 
edge radius. 


Figure 25 shows superimposed plots of alternate developed stream- 
surface sections. Figure 26 is a similar plot of all of the manufac- 
turing sections. 


ec. Location cf Stack Axis 


The stator stack axis is located at an axial coordinate of 
4,7250 inches, measured from the same origin as was used to define 
the annulus geometry. 


G. Root Fillet 
Between points 3/4 inch in from the leading edge and 1/2 inch 
in from the trailing edge, the fillet is 1/4 inch radius on both sides 


of the biade at both hub and casing. The fillet is smoothiy de- 
creased to 1/16 inch radius at the blade edges. 
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POINT VO XP yp XS. YS 

1 =24 33285 0273382 25336395 0275931 

2 ~2e30979 9475157 ~2,31571 0071564 

3 2428609 9672012 —22 239323 0.67283 

4 2626150 9655941 2226942 Oo 63267 

5 —2523537 0.55933 2924419 0. 59483 

5 2. c08iL 362996 2.21740 0.55740 

7 -2517936 0650110 -23 18891 0.5211% 

3 ~20e1489% 0557275 —-2215357 004260i 

) 2611667 1334487 Zs 12624 0045195 

19 —Z250&240 Os51742 -2,09177 0641290 

i il 2604595 02439037 -23 05569 Oo 38683 
2 23007146 0246369 -2¢01589 0.35579) ; 
| 13 1.96553 Qo 43733 ~i297408 0232550 3 
\ 14 1492229 0041144 -1.929383 0629624 E 
} 15 1.87637 0033594 1.89335 0625795 4 
a) 1.82831 0236092 1283457 0624067 3 
17 -1077827 0.33643 is 73409 0621444 ; 
13 1572642 0631253 —1, 73153 0613929 4 
| i3 1567294 9625927 -12 07745 0023527 qq 
yaa) ~1261303 023670 ~1,62198 0024240 7 2 

él 1056195 0324437 —1,56531 0.12071 

22 -1,.50441 0322334 le 50764 0219023 

: 23: 1246532 0 20363 -1544913 0608093 

{ 24 1.33818 0+13430 -12390i1 0005233 

i 25 Le 22906 0.16585 ~i» 33061 0.04023 

1 26 le 269% 0514337 1.27035 0.03076 

| 27 “1221011 = 9013183 -1,22105 0001655 
{ 20 “1515054 Os21l62T 1.215133 0.0036! E . 
i 29 2609135 0e10171 2309185 -0.00803 a 
39 1003237 0,038E5 -i203279  -0202851 E 
31 0297379 0207553 -0¢97399 0.02772 : 
; 32 -0.91568 0299401 -0,91575 -0.03573 * 
' 33 0s 35809 0405343 “U,E538L2 0004254 q 
| 34 -0450107 0204331 -0,80104 =0.04819 5 
35 0274402 3003516 -0:74456 0.05271 4 
| 35 0068876 0302744 —-0.605863 —-0.65619 a 
37 00632 +3 0302064 -0263339 0205840 : 
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win tterative Llirs reerticatin gprreedure was employed with the 
a lie, alffurkrn accc oitel rut wits lesz cetimistice assuxptiens for 
2 ..fé oilses im orfer t- estimate tne pretatile perfermance of the 
vy .rpllsa ctage anen Sirst tested. Tne 2cefficient cf the Mack numter 
*orctn £7 13) «ao imozreszced frer 2.4467 ta 1.1, which is equivalent 
%: trcuritg¢ t21a2 2 ncrctal encla steurs in tne rotor et the relative 
fret Vac. tatter and in tne itator at the exranded sucticn surface 
“st. G.faap. Tnaese walt iaticns may still te semechat cptimictic te 
* < setent that an elevated retror tic aiffusicn less waz not introe- 
g.le2. lay nalrulaticonz: were mage u:inse tne crogram cptisn of speci- 
fisarctir relative exit fli< angeles. Ecter deviaticn angles were 
B204T32 2 zeal t: their derign vaiuec. 

Tne pecults salruisteg at tne gecign flew rate cf 3).5 iscec 
ait: tae atove dercericed ohange ares 


Fctor trtal pressure ratio 3.342 
Eeter isertrepie efficiency 0. 2@9 
Stace total pressure ratic 2.738 
Stare isentrcpic efficiency 0.772 


“eridional velocity at the rclor exit on tre tip streamline dropped 
from 613.3 ft/see to $19.4 ft/eee and boundary layer Dicckage in the 
rotor exit plane rese from 0.0208 to §.0268 of the cress-sectional 
area. Tre relative total pressure less ccefficients resulting fren 
this analysis are shown for the rotor and stator in Fig 27. 


2. PREDICTED INTRA-SLADE-ROW PERFORMANCE 


The distributions of loss cbtained as described above and pre- 
sented in Fig 27 were used for an intra-blade-row analysis in the 
identical manner previously described in Section IV.2 for the criginal 
design. There were tyc ebjectives to this calculation. The first 
objective was to determine the choking fiow for the stage at loss 
levels more likely to occur in an initial test. The second objective 
was to determine to what extent the internal fiow distribution might 
change under these circumstances. 


With respect to the first objective, it was determined that the 


age would not pass the design fiow of 30.0 1b/sec but that it wouid 
$3 slightiy more than 29.0 lo/sec. Thus, 2 reduction of flow on the 
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orger cf terse percent Lelsu desien eam te anticinzated. Choxing fs 
treaicted te cccur at tre first caleulaticn station deenstrear of the 
~tator leading ¢lve. Tric fs in mntrast to tre Mach nurter distri- 
Guticn at decten ccs levels far wiien tre potcr, ratner than tre 
Stater, ecentrolc ezccine. Rotcr incidence angicc esrressonding to 
tne reduced fiew are arorecximately ane decree s.igter than design. 
Stator incidence angler are rredicted to te within one degree of de- 
Sigv: Values On @very ctreamcurface except tne cuter easing, where the 
rredicted value caeerpe up tc ten degrees atove desien. 

Concerning the cecend octeative, eines tne axial distribution af 
Static rressure was ehccen a: the carareter for enicn tne desion 
shouid te ontiniced, che statie créeccure dfetricuticn corresponding 
to the higher lose levels cnculd give an indleation of the departure 
from ideal ecenditicns. This ts cresented in Fig; c@, Corparing this 
with Fig 17 whien presents tne design distrizuticn, tne enceaed con- 

Stion cf the stator arrears net tc slerificantiy affect the stator 
@xiai pressure gradient cut if can te exrected ta cause tie rotor ts 
operate at a more throttled condition tran its reax ¢fficiency sper- 
ating point. Tnis would tyrically result in @ steerer axial cressure 
gradient in the entrance réricn of tre rotor fallesed ty a plateau 
as shown in Fir 28. Until steps are takzen tc reduces rotor icsses or 
to increase stator flicw area, it is unlinxels that the reoter will re 
aiie to eperate at its roint cf sean efficteney at desizn sreed. This 
eonelusicn is grediecatea upon the assurrticn trat maximum rotor ¢ffi- 
clency is likely tre oecur at or very near tne throttle coint at «nich 
the rotor just becomes unchosed. This operating condition tcrréesponds 
to the minimum axial rcrecsure gradient which can ie chtrained at a 
throttle point near the maximum rotor rressuré ratic. 
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